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SUMMARY 

This  is  a  preliminary  analysis  of  digitized  satellite  cloud  brightness 
data  that  will  be  used  in  developing  a  model  for  large-scale  analysis  over 
data-void  tropical  areas,  which  is  based  on  Wallace's  (1971)  postulations 
that  such  data  may  be  used  to  estimate  vertical  motion  in  the  tropics. 

The  relationship  between  satellite  cloud  brightness  and  large-scale 
200-mb  divergence  on  5  x  5  grids  in  the  tropical  western  North  Pacific  is 
examined  from  several  angles.   Correlation  coefficients  between  the  two 
fields  are  computed  on  a  daily  basis  for  the  region  of  study  and  at  each 
grid  point  for  8-month  time  series.   The  coefficients  are  mostly  low  although 
a  positive  correlation  is  generally  indicated.  The  low  correlation  may  be 
due  to  the  quality  of  the  available  data,  especially  for  the  divergence 
which  is  kinematically  computed  from  analyzed  winds.   Spectrum  analysis  is 
also  performed  to  determine  dominant  synoptic  time  scales  for  both  parameters 
in  the  period  of  study.  Two  common  period  bands,  one  centered  <~  10  days  and 
the  other  ~  5  days,  are  found.   Cross-spectra  between  the  two  parameters 
indicate  that  their  phase  differences  are  small  (<  1/4  cycle)  for  both  bands, 
but  the  coherence  squares  are  somewhat  lower  than  those  found  previously  by 
Wallace  for  the  4-  to  5-day  waves  using  direct  radiosonde  data.  Horizontal 
structures  of  both  parameters,  determined  from  inter-longitude  cross-spectra 
in  the  two  bands,  suggest  that  the  brightness  is  much  more  organized  than 
divergence,  but  the  two  fields  resemble  each  other  whenever  an  organized 
pattern  in  divergence  can  be  found. 

The  results  of  this  study  lend  some  support  to  Wallace's  postulation, 
but  also  indicate  some  problems  for  using  the  present  data.   It  is  antici- 
pated that  the  correspondence  with  the  actual  divergence  field  will  be 
greatly  improved  if  satellite  infrared  data  is  incorporated  into  the  bright- 
ness data. 
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1 .   Introduction 

Satellite  photographs  of  tropical  regions  have  revealed  that  much 
of  the  tropical  cloudiness  is  in  the  form  of  large,  connected  masses 
of  very  bright  clouds.  These  masses  are  composed  of  many  individual 
cumulus  clouds,  whose  extensive  cirrus  canopies  with  a  horizontal  scale 
of  several  hundred  km  or  larger  account  for  most  of  their  brightness. 
These  large,  bright  cloud  masses,  or  "cloud  clusters"  as  they  are  now 
called,  have  received  considerable  attention  in  recent  years.   Chang 
(1970)  and  Wallace  (1970,  1971),  using  time- longitude  series  of  satel- 
lite photographs,  found  that  many  cloud  clusters  over  the  tropical 
Pacific  are  well  organized  and  have  characteristics  resembling  those 
of  synoptic-scale  westward  propagating  wave  disturbances.   Further  evi- 
dence by  Reed  and  Recker  (1971)  showed  that  these  propagating  cloud 
clusters  in  the  western  North  Pacific  are  linked  to  synoptic-scale  waves 
in  the  wind  field.  Williams  and  Gray  (1973)  and  Yanai  et  al.  (1973) 
have  also  deduced  the  mean  properties  of  these  western  North  Pacific 
cloud  clusters  and  found  several  vertical  characteristics  of  the  active 
clusters  similar  to  those  found  by  Reed  and  Recker  (1971)  and  Nitta 
(1972)  for  the  wave  troughs. 

Due  to  the  poor  radiosonde  coverage  over  the  vast  tropical  oceans, 
it  is  natural  to  try  to  utilize  the  available  satellite  data  as  an  aid 
for  tropical  analysis.   Digitized  cloud  brightness  data,  acquired  by 
the  Advanced  Vidicon  Camera  Systems  of  the  operational  satellite  system, 
have  been  available  since  1967.   Because  of  the  relationship  between 
the  bright  cloud  clusters  and  synoptic-scale  tropical  motions,  this 
data  becomes  a  potential  source  for  such  a  purpose.   In  fact,  evidences 


of  relations  between  the  satellite  brightness  data  and  other  meteoro- 
logical parameters  have  been  found  by  several  investigators.   The 
correlation  between  brightness  and  precipitation  has  led  many  attempts 
to  estimate  rainfall  using  satellite  brightness  [for  example,  see  the 
review  by  Martin  and  Scherer  (1973)].   Wallace  (1971)  also  found  indi- 
cation of  good  correspondence  between  area-averaged  cloud  brightness 
and  synoptic-scale  vertical  motion  associated  with  the  4-  to  5-day 
tropical  disturbances  in  the  Kwajalein-Eniwetok-Ponape  triangle  of  the 
Marshall  Islands  and  the  Guam-Truk-Yap  triangle  farther  west.   He 
postulated  that  it  may  be  possible  to  obtain  an  estimate  of  the  verti- 
cal motion  field  over  the  tropics  from  the  cloud  brightness  data  alone. 
If  this  turns  out  to  be  the  case,  the  satellite  cloud  brightness  data 
will  certainly  be  one  of  the  best  available  indicators  of  disturbed 
weather  in  the  tropics. 

There  are  physical  reasons  to  think  that  brightness  data  is  closely 
related  to  the  large-scale  vertical  motion  field.  The  bright  cirrus 
canopies  associated  with  the  cloud  clusters  "are  produced  by  outflow 
from  and  remnants  of  the  cumulonimbi.   In  general,  developing  and  con- 
servative clusters  maintain  their  cumulonimbi  from  a  steady  low-level 
mass  convergence;  clusters  gradually  die  when  their  low-level  mass 
convergence  is  eliminated"  (Williams  and  Gray,  1973).   Because  of  the 
readily  available  low-level  moisture,  the  large-scale  upward  motion 
provides  a  favorable  environment  for  the  development  and  existence  of 
the  clusters.   Furthermore,  brighter  clouds  are  usually  thicker  and 
therefore  may  indicate  larger  amount  of  latent  heat  release  (Gruber, 
1974).   Because  the  magnitude  of  the  latent  heating  is  nearly  one  order 
greater  than  the  temperature  fluctuations  at  most  levels  for  large-scale, 


convectively-active  tropical  motions,   thermal  energy  balance  requires 
that  the  large-scale  vertical  motion  must  be  approximately  proportional 
to  heating  (Wallace,  1971;  Holton,  1972).  Thus  one  could  expect  a  posi- 
tive correlation  between  cloud  brightness  and  large-scale  vertical  motion 
by  reason  of  heating  alone. 

The  correspondence  between  brightness  and  vertical  motion  may  also 
be  extended  to  large-scale  horizontal  divergence  at  levels  where  diver- 
gence is  large,  since  many  observational  studies  [Wallace  (1971),  Reed 
and  Recker  (1971),  Nitta  (1972),  Williams  and  Gray  (1973),  Yanai  et  al. 
(1973),  Reed  and  Johnson  (1974)]  have  produced  very  similar  profiles  of 
the  vertical  motion  (or  divergence)  over  tropical  areas  where  convection 
is  active.   The  same  conclusion  can  also  be  inferred  from  the  various 
vertical  heating  profiles  shown  by  Hayashi  (1974)  from  analyses  of  tropi- 
cal output  from  the  general  circulation  model  of  the  NOAA  Geophysical 
Fluid  Dynamics  Laboratory  at  Princeton.   In  addition,  Williams  and  Gray 
(1973)  have  shown  a  clear  in-phase  relation  between  cloud  clusters  and 
upper- level  divergence  in  their  composite  study. 

Due  to  the  constraint  of  available  radiosonde  network,  the  study  of 
correlation  between  brightness  and  vertical  motion  by  Wallace  (1971)  is 
limited  to  just  two  "grid  points"  --the  Kwajalein-Eniwetok-Ponape  and 
the  Guam-Truk-Yap  triangles.   It  will  be  quite  significant  if  this  posi- 
tive correlation  between  brightness  and  vertical  motion  or  upper-level 
divergence  can  be  established  over  a  broader  area.   For  example,  since 
deep  tropical  motions  are  usually  quasi-barotropic,  it  may  be  possible 
to  represent  the  divergence  by  the  brightness  data  in  a  one-level  vorti- 
city  equation  to  obtain  some  kind  of  solution  which  may  resemble  the 
large-scale  flow  field.   This  can  then  be  used  as  a  first  estimate  for 


diagnostic  purposes  over  regions  where  wind  observations  are  sparse  or 
void.   The  possibility  of  this  kind  of  application  has  been  increased 
recently  by  the  available  satellite  infrared  data  which,  by  supplying 
information  of  cloud-top  temperatures  and  therefore  cloud-top  heights, 
may  be  used  to  modify  the  brightness  data  and  improve  the  estimate  of 
vertical  motion  and  upper-level  divergence  fields. 

A  model  of  using  satellite  brightness  data  based  on  the  foregoing 
discussion  is  currently  being  developed  at  the  Naval  Postgraduate  School. 
This  study  is  a  preliminary  analysis  of  satellite  brightness  data  which 
will  be  used  in  that  model.  Our  primary  purpose  is  to  examine  the 
degree  of  correspondence  between  the  brightness  data  and  the  large-scale 
vertical  motion-divergence  field,  using  the  analyzed-wind  data  from  the 
National  Meteorological  Center.  This  data  is  available  at  all  tropical 
grid  points  and  thus  enables  us  to  examine  a  large  area  of  tropical 
western  Pacific.   In  addition  to  the  4-  to  5-day  period  disturbances 
studied  by  Wallace,  we  will  examine  the  relationship  between  the  two 
fields  from  several  angles:  as  functions  of  space,  time,  and  dominant 
frequency  bands.   Correlation  and  spectrum  analyses  will  be  used.  The 
plan  of  developing  the  model  for  using  satellite  brightness  data  as  an 
aid  to  large-scale  analysis  over  data-void  tropical  regions  will  also 
be  briefly  discussed  in  the  final  section. 
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2.   Data  sources  and  periods  of  study 

The  general  area  selected  for  this  study  is  the  tropical  western 
North  Pacific,  from  the  equator  to  25N  and  125E  to  175E.  There  are  two 
reasons  for  this  selection:   Firstly,  previous  spectral  wave  studies 
and  satellite  studies  have  indicated  that  this  is  an  area  of  strong 
cloud  cluster  activity.   Secondly,  this  region  has  eleven  rawinsonde 
stations  including  the  Caroline  and  Marshall  Islands,  which  makes  the 
analyzed  wind  data  relatively  more  reliable  compared  to  other  areas  in 
the  tropics.   Fig.  1  shows  the  region  and  the  rawinsonde  network. 

The  digitized  cloud  brightness  and  upper-level  wind  data  used  in 
this  study  were  provided  by  the  National  Center  for  Atmospheric  Research 
(NCAR) .   The  cloud  brightness  data,  acquired  by  NCAR  from  the  National 
Environmental  Satellite  Service  (NESS),  are  digitized  values  averaged 
for  5x5   latitude-longitude  squares.   Upper-level  component  winds  at 
the  700,  500,  300,  250  and  200  mb  levels  were  acquired  by  NCAR  from  the 
National  Meteorological  Center  tropical  grid  analyses.  These  winds  are 
given  at  grid  points  approximately  coinciding  with  those  for  the 
brightness  data,  and  are  the  result  of  analyses  based  on  all  available 
rawinsondes,  aircraft  observations  and  a  few  winds  deduced  from  cloud 
drift  as  seen  from  satellite  photographs.  NCAR  also  provided  rawinsonde 
observations  for  the  region  of  study,  which  were  used  to  obtain  analyzed 
wind  fields  at  low  levels.   In  all  cases,  the  0000  GMT  wind  data  are 
considered  as  being  closest  to  the  time  of  satellite  observation. 

Two  periods  have  been  selected  for  this  study:   July  1969  and  19  April 
-  20  December  1971.   The  first  period,  which  is  just  a  randomly  selected 
summer  month,  was  first  used  to  compute  daily  correlations  between  cloud 
brightness  and  vertical  motion  fields.  The  second  period  was  selected 
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when  we  decided  to  extend  the  period  of  analysis  and  to  resolve  corre- 
lations for  dominant  frequency  bands.   By  that  time  we  had  learned  of 
the  availability  of  the  NOAA-1  satellite  infrared  data  digitized  by  NESS, 
which  covers  approximately  a  three-month  period  starting  19  April  1971. 
We  hope  that  this  infrared  data  can  be  used  at  a  later  date  to  modify 
and  compare  with  the  present  results. 
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3.   Results 

3.1  Daily  correlations  for  July  1969 

The  period  of  July  1969  contains  23  days  of  useful  brightness  data, 
which  were  used  to  compare  with  the  divergence  and  vertical  motion  (pres- 
sure velocity  od)  on  a  daily  basis.  The  latter  fields  were  computed 
kinematically  on  the  5x5  grids  subjected  to  the  boundary  constraints, 
oo  =  0  at  1000  mb  and  100  mb.  The  correction  necessary  to  satisfy  these 
constraints  is  distributed  linearly  with  respect  to  pressure,  resulting 
in  an  adjusted  divergence  and  vertical  velocity  at  each  level.  The 
validity  of  such  calculations  on  a  5  x  5  grid  is  undoubtedly  question- 
able, and  even  more  so  when  one  considers  that  our  wind  data  is  already 
smoothed  by  the  analysis  procedure.  However,  thus  far  this  is  the  only 
data  available  over  a  broad  tropical  oceanic  area;  and,  for  our  purpose, 
we  believe  that  a  crude  representation  of  the  vertical  motion  or  diver- 
gence field  will  be  adequate  to  show  at  least  whether  the  correlations 
with  brightness  found  by  Wallace  (1971)  on  a  time  series  basis  at  two 
"grid  points"  can  be  extended  to  a  larger  area  on  a  daily  basis.  Further- 
more, the  computed  divergence  profiles  in  areas  of  high  brightness  usually 
show  a  deep  layer  of  convergence  capped  with  a  relatively  shallow  layer 
of  divergence,  which  are  generally  in  good  agreement  with  the  results 
found  by  Wallace  (1971),  Reed  and  Recker  (1971),  Williams  and  Gray  (1973) 
and  others,  thus  lending  some  confidence  to  the  calculations. 

The  daily  correlation  coefficients  of  brightness  vs  vertical  velocity 
and  divergence  were  computed  at  all  levels.   Generally,  the  300-  to  400-mb 
vertical  velocities  and  the  200-mb  divergence  show  better  correlations 
with  brightness  as  compared  t&   other  levels.  This  is  somewhat  expected 
because  previous  observational  studies  all  found  that  the  non-divergent 
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level  is  ~  300-  to  400-mb,  while  the  200  mb  level  is  the  outflow  level 
with  maximum  divergence.   Table  1  lists  the  daily  correlation  coeffi- 
cients for  the  300-mb  vertical  velocity  (^30o^  and  200"mb  divergence 
(div-,-^) .   The  sample  size  is  66,  which  gives  us  a  value  of  0.24  for  the 


95$  significant  level  and  0.20  for  the  90%  significant  level,  using 
Panofsky  and  Brier's  (1965)  formula.  Out  of  the  23  days  of  useful  data, 
14  days  show  a  significant  correlation  at  90$  level  or  better  between 
brightness  and  ^^q^   »  with  13  days  negative  and  1  day  positive.  These 
correlation  coefficients  are  quite  low  as  compared  to  that  (inferred 
from  coherence  square)  found  by  Wallace  (1971)  for  the  4-  to  5-day  waves 
at  the  Rwajalein-Eniwetok-Ponape  triangle,  although  they  do  indicate  a 
tendency  of  in-phase  relationship  between  the  two  parameters  on  a  daily 
basis.   It  is  possible  that  the  lower  correlation  is  partly  due  to  the 
quality  of  our  computed  divergence  data  which  is  based  on  analyzed  winds, 
but  our  results  do  suggest  that  the  correspondence  between  the  two  fields 
is  far  from  perfect.   In  column  2  of  Table  1  the  correlation  coefficient 
between  ^onn     and  a  smoothed  brightness  data  is  also  listed.  Here 
the  smoothing  procedure  is  done  by  a  simple  5-point,  centrally-weighted 
spatial  averaging  in  order  to  reduce  the  effect  of  isolated  high  or  low 
spots  that  appear  occasionally.  However,  this  averaging  turned  out  to 
be  of  little  use.  Although  most  of  the  correlations  are  slightly  higher, 
they  are  not  too  different  from  those  for  the  original  brightness. 

The  correlation  of  brightness  vs  div„on  ,  as  shown  in  Table  1,  seems 
quite  similar  to  that  for  <^30q  although  many  coefficients  are  slightly 
lower.  This  similarity  in  correlation  with  brightness  is  somewhat 
expected  from  the  previous  observation  that  the  vertical  profile  of 
vertical  velocity  (or  divergence)  remains  approximately  the  same  in  the 
regions  of  strong  cloud-cluster  activities. 
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In  addition  to  the  correlation  calculations,  we  also  visually  compared 
the  plotted  contour  patterns  of  the  brightness  with  those  for  the  co_nn 
and  div„00  fields.  This  subjective  examination  reveals  that  the 
area  enclosed  by  10  N,  20  N,  130  E  and  155  E  generally  shows  better  cor- 
respondences between  the  different  parameters.   In  fact,  the  eastern 
portion  of  the  region  (160  E  -  175  E)  possesses  comparatively  few  organized 
brightness  patterns  in  the  contoured  diagrams  and  usually  exhibits  a 
widespread  area  of  zero  or  low  brightness.  This  increase  in  organization 
and  intensity  of  brightness  toward  the  western  portion  may  be  a  reflection 
of  the  westward  intensification  of  wave  disturbances  in  this  region  found 
by  Chang  et  al.  (1970)  and  Reed  and  Recker  (1971).   Chang  (1970)  and 
Sikdar  et  al.  (1972)  also  noted  a  marked  geographical  variation  in  cloud 
activity  over  the  Pacific. 

3.2  Correlations  for  1971 

Because  the  July  1969  correlation  results  are  similar  for  o)-nn  and 
div_nn  ,  we  decided  to  use  the  horizontal  divergence  only  for  correla- 
tion with  the  cloud  brightness  for  1971.  The  divergence  is  derived 
directly  from  the  horizontal  wind  components  at  all  levels  without  any 
adjustment  by  vertical  boundary  constraints.  The  200-mb  level  again  turns 
out  to  be  the  maximum  divergence  level  for  almost  all  the  days  in  this 
period,  so  only  results  using  div_n„  are  presented  here  and  in  3.3. 

Table  2  lists  the  daily  correlation  coefficients  of  brightness  vs 
div_nn  for  a  four-month  summer-autumn  period,  namely  July-October  1971. 
Among  the  119  days  with  useful  data,  45  days  show  positive  correlation 
at  the  90$  level  or  better,  while  4  days  show  negative  correlation. 
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These  results  are  similar  to  those  of  July  1969  and  also  provide  a  weak  indi- 
cation of  positive  correlation  between  the  two  fields  on  a  daily  basis. 

The  correlation  coefficients  between  the  8-month  (19  April  -  20  December 
1971)  time  series  of  each  parameter  were  also  computed  at  all  grid  points 
in  order  to  examine  the  spatial  distribution  of  the  correlations.  The  results 
are  listed  in  Table  3.  Among  the  66  grid  points  for  the  entire  region,  25 
points  show  positive  correlation  at  the  90$  level  or  better  vs  3  points  with 
negative  correlation.  This  ratio  is  compatible  with  the  daily  computations. 
In  fact,  all  three  negative  points  are  in  the  southeast  corner  of  the  region, 
and  most  interior  points  show  positive  correlations.  This  situation  is  some- 
what similar  to  the  results  of  our  visual  inspection  of  the  July  1969  data 
mentioned  earlier.   It  thus  seems  to  us  that  if  a  region  bounded  by  5N,  20N, 
130E  and  160E  is  chosen,  the  daily  correlations  would  be  much  improved. 

3.3  Spectrum  analysis  for  1971 

The  correlation  coefficients  between  brightness  and  div9f)f)  listed 
in  Tables  1-3  are  quite  low  as  compared  with  the  coherence  for  the  4- 
to  5-day  waves  found  by  Wallace  (1971)  in  the  Kwajalein-Eniwetok-Ponape 
triangle.  Therefore,  it  may  be  worthwhile  to  do  a  spectrum  analysis 
of  our  data  and  compare  the  results  with  Wallace's.   In  addition,  the 
existence  and  importance  of  synoptic-scale  disturbances  in  the  tropics 
have  been  well  recognized.   It  is  reasonable  to  expect 


However,  because  of  the  many  uncertainties  in  the  computed  divergence 
data,  the  percentage  of  positive  correlation  in  Tables  1-2  should  not  be 
viewed  as  equivalent  to  the  probability  of  finding  a  positive  correlation 
between  brightness  and  the  actual  div„nn  on  a  given  day. 
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that  the  variations  in  both  the  divergence  and  the  brightness  fields 
are  highly  influenced  by  these  disturbances,  and  these  variations  may 
be  characterized  by  certain  dominant  time  scales  associated  with  the 
synoptic-scale  disturbances  (which  may  or  may  not  be  related  to  the 
easterly  waves).  However,  it  is  unclear  whether  our  previous  correla- 
tion results  are  mostly  due  to  the  very  large  scale  controls  (such  as 
monsoonal  circulations,  location  of  ITCZ,  etc.)  or  the  synoptic  dis- 
turbances.  If  we  want  to  try  to  use  the  brightness  data  to  represent 
divergence  or  vertical  motion  on  a  daily  basis,  it  is  necessary  that 
the  characteristic  behavior  of  the  two  parameters  resemble  each  other 
in  these  synoptic  time  scales.   For  these  reasons,  spectrum  and  cross- 
spectrum  analyses  are  performed  for  the  19  April  -  20  December  1974 
time  series  of  both  parameters  at  all  66  grid  points  within  the  region 
of  study. 

Before  the  spectral  results  are  presented,  it  is  useful  to  briefly 
discuss  our  aim  as  compared  to  those  of  other  spectrum  studies  of 
satellite  brightness  data  attempted  previously  by  tropical  meteorolo- 
gists.  There  have  been  a  number  of  such  studies  (Tanaka  and  Ryuguji, 
1971,  1973;  Murakami  and  Ho,  1972;  Sikdar  et  al . ,  1972,  and  others) 
carried  out,  primarily  for  the  purpose  of  identifying  tropical  wave 
disturbances.  Although  some  have  yielded  interesting  results,  Wallace 
and  L.  Chang  (1972)  have  suggested  that  the  number  of  spectral  peaks 
found  by  these  studies  may  be  so  numerous  that  it  appears  rather  dif- 
ficult to  extract  useful  information  from  the  results.  They  have 
also  pointed  out  that  the  time  spectra  of  the  present  brightness  data 
are  highly  variable  in  both  space  and  time,  and  even  neighboring  grid 
points  may  exhibit  quite  different  spectral  characteristics.  However, 
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our  primary  purpose  in  using  spectrum  analysis  is  not  to  identify  general 
wave  scales,  but  rather  to  determine  how  much  the  brightness  data 
resembles  the  divergence  when  both  parameters  are  decomposed  into 
various  time  scales  which  are  of  synoptic  interst  in  our  period  of 
study. 

All  the  8-month  time  series  are  reduced  to  a  length  of  198  days 
after  the  application  of  a  high-pass  filter,  which  does  not  have  much 
effect  for  periods  <  20  days.   (The  response  at  the  20-day  period  is 
97$.)  A  lag  of  20  days  is  used  for  all  computations.   Figs.  2  and  3 
are  the  power  spectra  at  each  grid  point  for  brightness  and  div_00  , 
respectively.  After  examining  these  spectra,  we  found  that  two 
frequency  ranges,  .075  -  .125  cycle  day  '  (or  period  8-13.3  days)  and 
.175  -  «125  cycle  day   (or  4.4-5.7  days),  contain  substantial  portions 
of  the  total  variance  of  most  series,  hence  it  is  convenient  to  choose 
these  two  bands  for  comparing  the  behavior  of  the  two  parameters. 
For  simplicity,  they  shall  be  subsequently  called  the  10-  and  5-day 
bands,  respectively.  Although  we  are  not  especially  concerned  about 
the  location  of  the  spectral  peaks,  it  is  nevertheless  noteworthy  that 
most  of  the  major  peaks  do  fall  into  these  two  bands,  and  that  previous 
studies  by  Sikdar  et  al.  (1972),  Tanaka  and  Ryuguji  (1973)  and  Gruber 
(1974)  have  all  found  similar  two-band  structures  in  the  spectra  of 
tropical  brightness  data.   It  thus  seems  to  us  that  although  there  is 
a  high  variability  in  the  brightness  spectra,  simple  interpretation 
of  a  two-band  distribution  appears  to  be  quite  consistent  in  several 
studies  and  should  not  be  overlooked.  We  also  notice  that  a  3-day 
spectral  peak  appears  occasionally, but  they  seldom  appear  concurrently 
in  both  parameters  and  the  fraction  of  total  variance  involved  is 
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relatively  small  in  most  cases. 

The  coherence  square  and  phase  difference  between  brightness  and 
div_00  at  each  grid  point,  averaged  within  each  dominant  frequency 
band,  are  listed  in  Table  4.  The  averaging  process  results  in  a  degree 
of  freedom  of  ~  25,  which  gives  a  coherence  square  of  0.17  at  the  99f0 
significant  level,  0.12  at  the  95$  level  and  0.09  at  the  90$  level. 
Only  those  values  significant  at  the  90$  level  or  better  are  shown 
in  Table  4.   For  the  10-day  band,  24  grid  points  have  significant 
coherence  squares,  and  most  of  them  show  phase  differences  from  in- 
phase  to  ~  t   cycle  between  the  two  parameters.   For  the  5-day  band, 
20  points  have  significant  coherence  squares.  With  the  exception  of 
3  points  (5N-135E,  20N-145E  and  20N-150E),  the  phase  differences  also 
range  between  in-phase  and  —  cycle  between  the  two  parameters,  but 
there  appears  to  be  more  t   cycle  phase  differences  in  this  band  than 
in  the  10-day  band.   Comparison  of  our  results  for  both  frequency 
bands  with  those  of  Wallace  (1971)  for  the  4-  to  5-day  waves  showed 
that,  although  the  two  parameters  are  much  more  in-phase  than  out-of- 
phase,  the  coherence  squares  were  lower  and  the  phase  differences 
are  larger.  This  is  perhaps  again  due  to  our  use  of  smoothed, 
analyzed  winds  for  computing  the  divergence  field. 

Another  indication  of  correspondence  between  the  two  parameters 
in  the  5-  and  10-day  bands  may  be  found  by  comparing  the  horizontal 
structures  deduced  from  inter- longitude  cross-spectra  for  both  para- 
meters at  each  latitude.  The  results  are  shown  in  Figs.  4-5.   For 
each  latitude  the  150E  series  is  used  as  the  base  series  to  cross 
with  other  series  at  the  same  latitude.  Each  phase  difference  is 
plotted  against  the  longitudinal  distance  from  the  base  series,  and, 
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unless  the  points  are  widely  scattered,  a  best  fitting  straight  line 
is  drawn  for  each  latitude  weighted  by  the  significant  level  of  the 
coherence  square  associated  with  each  phase  value.   This  technique  has 
been  commonly  used  in  spectral  wave  studies  for  determining  horizontal 
wavelengths  and  propagating  directions. 

For  the  10-day  band,  a  westward  movement  in  the  brightness  data 
with  a  "wavelength"  of  ~  50  can  be  inferred  at  all  latitudes  except 
the  equator,  although  this  simple  interpretation  seems  to  be  inappli- 
cable in  the  eastern  portion  of  ION  and  25N,  and  the  easternmost 
section  of  20N.  The  divergence  field  is  less  organized,  but  a  similar 
structure  does  appear  at  5N,  15N  and  20N,  the  three  latitudes  where 
the  brightness  results  are  best.   In  addition,  another  similarity  can 
be  found  at  ION.  The  divergence  series  are  in-phase  between  140E-175E 
and  the  brightness  series  are  also  in-phase  between  145E-170E.  Both 
results  indicate  an  east-west  orientation  and  may  suggest  the  presence 
of  an  ITCZ  with  its  intensity  fluctuating  around  a  period  of  r^  10  days. 

The  structures  for  the  5-day  band  seem  to  be  less  regular  as  com- 
pared to  the  10-day  band.  Here  the  brightness  is  again  more  organized 
than  the  divergence.  At  least  a  partial  best  fitting  line  can  be 
drawn  for  brightness  at  all  latitudes,  although  the  "wavelength"  range 
from  ~  25  at  the  northern  latitudes  to  ~  50  at  the  southern  latitudes. 
Only  two  organized  patterns,  at  5N  and  20N,can  be  found  for  the 
divergence,  but  in  each  case  they  do  remarkably  resemble  those  of  their 
respective  brightness  fields. 

The  fact  that  the  brightness  data  is  more  organized  in  Figs.  4-5 
reaffirms  our  belief  that  its  quality  is  better  than  that  of  the  diver- 
gence, which  is  expected  as  mentioned  before.   Despite  the  noise  present 
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in  our  divergence  field,  whenever  an  organized  pattern  in  the  divergence  can 
be  found,  it  resembles  the  brightness  pattern.   This  may  be  viewed  as  another 
indication  of  the  correspondence  between  the  two  fields. 
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4.   Concluding  remarks 

In  this  study  we  have  examined,  from  several  angles,  the  relation- 
ship between  the  satellite  cloud  brightness  and  the  divergence-vertical 
motion  field  on  the  5x5  grid  in  the  tropical  western  North  Pacific. 
Our  results  appear  to  be  somewhat  mixed.   On  one  hand  the  correlation 
coefficients  and  coherence  squares  between  the  two  fields  are  consider- 
ably lower  than  those  previously  found  by  Wallace  (1971)  for  the  5-day 
waves  using  direct  radiosonde  data.   On  the  other  hand,  when  the  results 
of  various  comparisons  are  integrated,  they  do  form  a  consistent  picture 
of  positive  correlation  between  the  two  fields.   In  this  regard  the 
present  results  lend  some  support  to  Wallace's  contention  that  it  may 
be  possible  to  estimate  the  vertical  motion  or  upper-level  divergence 
field  over  the  tropics  from  digitized  satellite  brightness  data. 

Our  results  are  undoubtedly  influenced  by  the  quality  of  divergence 
data  which  are  deduced  from  analyzed  winds.  However,  until  the  obser- 
vational network  in  the  tropics  can  be  improved,  it  will  be  very 
difficult  to  obtain  better  divergence  values  from  wind  reports.   This 
is,  of  course,  the  precise  reason  for  the  attempts  to  use  satellite 
data.  There  are  also  problems  in  using  the  brightness  data  to  repre- 
sent cloud  clusters,  but  the  infrared  data  now  available  provides  a 
great  possibility  for  improvements.   Infrared  data  supplies  information 
on  cloud  top  temperatures  and  can  be  used  as  an  indicator  of  the  verti- 
cal extent  of  cumulus  development  and  therefore  the  vertical  extent 
of  latent  heat  release.   It  is  thus  possible  to  modify  the  brightness 
data  by  the  infrared  data  and  obtain  better  correlations  with  the 
divergence  field. 
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Moreover,  since  presently  there  is  no  reliable  direct  way  to  obtain 
divergence  or  vertical  motion  in  the  data-sparse  area  of  tropics, 
satellite  data  remains  as  one  of  the  best  indicators  of  vertical  motion. 
In  this  study  we  have  shown  that  it  is  possible  to  deduce  organized 
200-mb  divergence  patterns  associated  with  synoptic-scale  disturbances 
from  brightness  data.  There  may  also  be  other  applications  of  the 
satellite  data.  One  possible  use,  which  is  currently  being  investigated 
at  the  Naval  Postgraduate  School,  is  related  to  the  barotropic  and 
quasi-nondivergent  properties  of  the  large-scale  tropical  flow.  As 
mentioned  in  the  introduction,  it  may  be  feasible  to  describe  the  tropi- 
cal flow  field  at  a  single  level  by  a  barotropic  vorticity  equation, 
such  as  the  simple  diagnostic  model  by  Holton  and  Colton  (1972)  for 
the  200-mb  time  mean  circulations  during  the  northern  summer.  Using 
the  observed  time  mean  divergence  and  mean  zonal  wind,  they  success- 
fully reproduced  the  200-mb  mean  motion  field  by  the  steady  state  solu- 
tion of  a  linear  barotropic  vorticity  equation  with  a  simple  para- 
meterization of  the  strong  damping  process.   If  satellite  data  can  be 
used  to  represent  horizontal  divergence  at  200  mb,  a  similar  model  may 
also  be  developed  for  a  data-void  tropical  region  on  a  daily  basis. 
Here  the  mean  zonal  wind  may  be  represented  by  climatology  and  the 
boundary  conditions  determined  by  the  data-rich  surrounding  areas. 
Since  the  observational  results  by  Williams  and  Gray  (1973)  and  many 
others  have  shown  that  the  local  change  of  vorticity  is  rather  small 
compared  to  other  terms  in  the  vorticity  equation,  a  steady-state 
solution  obtained  from  such  a  simple  model  may,  to  some  degree,  resemble 
the  actual  flow  if  the  strong  damping  process  at  200  mb  (Reed  and 
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Recker,  1971;  Williams  and  Gray,  1973)  can  somehow  be  represented. 
If  such  a  model  proves  to  be  feasible,  it  will  be  very  useful  for 
obtaining  a  first  estimate  of  the  200  mb  flow  over  a  data-void  region 
in  the  tropics . 
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Table  1.   Daily  correlation  coefficients  in  hundredths  of 

satellite  cloud  brightness  vs  ^0™  anc*  div___  ^or 
July  1969.  Values  shown  in 

parentheses  are  significant  between  the  90$  and 
95$  levels;  values  shown  without  parentheses  are 
significant  at  the  95$  level  or  better. 

Smoothed 
Brightness   Brightness   Brightness  vs 

Date         VS  ^300     VS  ^300       div200 


1  July  1969       -  - 

2  - 

3  -24         -24  (22) 
6              (-21)        (-22)  (20) 

8  - 

9  ... 

10  ... 

11  - 

12  (20)         23  (-21) 

13  - 
15 
17 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 


- 

-25 

- 

-23 

-27 

24 

-23 

-23 

(21) 

-24 

-24 

26 

(-22) 

(-19) 

- 

-30 

-33 

25 

-32 

-41 

33 

(-22) 

-26 

(20) 

-32 

-30 

31 

-36 

-40 

36 

-60 

-62 

42 

-34 

-33 

27 

(-19) 

-30 

- 

25 


Table  2.  Daily  correlation  coefficients  in  hundredths  of 
satellite  cloud  brightness  vs  div9f)f)  for  July  - 
October  1971.   Values  shown  in 


parentheses  are  significant  between  the  90%  and 
95$  levels;  values  shown  without  parentheses  are 
significant  at  the  95$  level  or  better.  M  repre- 
sents missing  data. 


Day 

July 

August 

September 

October 

1 

- 

M 

- 

- 

2 

- 

28 

- 

- 

3 

(-22) 

(22) 

- 

- 

4 

- 

29 

- 

26 

5 
6 

7 

- 

28 

26 

M 

(22) 

_ 

31 

_ 

8 

28 

M 

42 

33 

9 

(22) 

-31 

(20) 

- 

10 

- 

- 

34 

36 

11 

32 

- 

- 

34 

12 

38 

M 

- 

(20) 

13 

30 

■» 

46 

M 

14 

34 

- 

- 

- 

15 

- 

- 

44 

- 

16 

- 

- 

52 

- 

17 

- 

- 

42 

- 

18 

- 

- 

(23) 

- 

19 

M 

- 

- 

- 

20 

- 

- 

36 

- 

21 

(23) 

(20) 

25 

- 

22 

40 

29 

- 

- 

23 

41 

- 
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24 
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(21) 

- 
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- 

25 
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(22) 

28 

26 
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(23) 

(21) 

36 

27 

- 

- 

- 
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28 

- 
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29 

29 
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(-22) 

- 

30 

39 

- 

- 

- 

31 
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i 

!! 

2V*     150£ 

1 

0- 

V 

m~r*   of  digitized  satellite  brightness  data  at 
Fig-  §.   Power  spectra  of  ^tize        variance  per  unit  frequency 
the  5°  x  5°2grid  ^%'^JT^T0-\   where  B  is  brightness  unit 
SSTU.  £■ ^  to'lO)  S  the  ahsclssa  is  period  (days). 
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Fig.    2.    (continued) 
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?       •  20  10    7     5    4         3 


7     5     4         J  2       -JO  10   7     5 


J      •  20  10    7    5 


Fig.     3 


Same  as   Fig.    2   except   for  200-mb  divergence  and  the  unit 
for  the  ordinate   is   sec"2  per  2tt  x  40"-1-  day"1  x   10"      . 
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Fig.    3.    (continued) 
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